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Abstract. The advent of the Spitzer Space Telescope has revolutionized 
our understanding of the history of star formation and galaxy mass assembly 
in the Universe. The tremendous leap in sensitivity from previous mid-to-far 
IR missions has allowed Spitzer to perform deeper, and wider, surveys than 
previously possible at these wavelengths. In this review I highlight some of the 
key results to come out of these surveys, and the implications these have for 
current models of galaxy formation and evolution. 



1. Introduction 

Spitzer, with its complementary suite of near-to-far infrared instruments, has 
made important contributions to virtually every aspect of astrophysics in the 6 
years since its launch. To satisfactorily review all of its accomplishments, even 
in the single field of galaxy formation and evolution, would be a tremendous 
task. Thus this all-too-brief review will cover only selected highlights of results 
coming from some of the key extragalactic legacy surveys, in particular; the 
Great Observatories Origins Deep Survey (GOODS; Dickinson et al. in prep), 
the Spitzer Wide-area InfraRed Extragalactic survey (SWIRE; Lonsdale et al. 
2003), The Far Infrared Deep Extragalactic Legacy survey (FIDEL; PI Dick- 
inson) and Spitzer observations taken as part of the Cosmic Evolution Survey 
(COSMOS; Scoville et al. 2007). In particular a focus on statistical measures, 
such as the luminosity function or two-point correlation function, is made as 
these are areas in which Spitzefs ability to survey large areas to exceptional 
depth has made the biggest impact. 

This review is broken down into three distinct parts; 1. the star forma- 
tion history of galaxies, 2. Properties of star forming galaxies, and 3. The 
environment of star forming galaxies, as seen by Spitzer. 



2. Star formation history of Galaxies 
2.1. Cosmic Infrared Background 

The most direct, albeit crude, measure of the mid-to-far IR properties of galaxies 
is to simply measure the total emission across the sky, which forms the Cosmic 
Infra-red Background (GIB; Hauser et al. 1998). While some of this emission 
comes from AGN, the bulk of this can be directly attributed to star forming 
processes in distant galaxies. Spitzer has enabled, for the first time, a large 
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fraction of the galaxies which make up the CIB at different wavelengths to be 
resolved. At 24/xm Papovich et al. (2004) directly observed the turnover in the 
number counts at ~ 0.2mJy and hence directly estimate a lower limit to the 
CIB at 24/xm to be i^U = "^-^to^ nW m"^ g^-i^ consistent with COBE DIRBE 
upper limit of vl^ < 7 nW m^^ sr^^ (Kashlinsky & Odenwald 2000; Hauser 
k Dwek 2001). Of this background ul^ = 1.9 ± 0.6 nW m'^ sr^^ (~70%) is 
contributed by sources > 60//Jy which are directly resolved by Spitzer MIPS. 

However even with the improved sensitivity of Spitzer a large fraction of 
the 70 &: 160 /im background cannot be resolved directly into individual sources. 
Dole et al. (2004) show that only ~20% (~10%) of the CIB at 70/xm (160/xm) is 
contributed by sources resolved by MIPS surveys at those wavelengths. One way 
to get around this is to "stack" the images at 70 (& 160) /im based on the known 
positions in shorter wavelength data. Dole et al (2006) performed stacking of the 
MIPS 70 & 160 /im images in three fields (CDFS, HDFN & Lockman Hole) based 
on the known positions of 24/im sources. The 24;um sources were distributed 
into 20 bins of descending flux density. Subsequently the 70 & 160 //m micron 
MIPS images were co-added at the positions of the 24;um sources in each flux 
density bin. In this way the contribution to the 70/Ltm and 160/im background as 
a function of 24//m flux can be assessed. Again using the observed behaviour of 
the 24^m number counts a robust lower limit on the 70/xm & 160/xm background 
of 7.1±1.0 & 13.4ibl.7 nW m'^ sr'^ Sources resolved at 24/xm (i.e. > 60/LtJy) 
contribute ^ 80% of this background at 70 Sz 160 /xm. 

2.2. Number counts &; Luminosity function of star forming galaxies 

One of the greatest surprises to come out of the earliest Spitzer observations 
was the great discrepancies between the observed number counts at MIPS wave- 
lengths and those predicted by prc-Spitzer IR population models. In particular 
Papovich et al. (2004) found that the peak in the 24^m counts was at a much 
lower flux density than predicted from phenomenological constrained to match 
ISO 15/xm surveys. Papovich et al., and subsequent studies, attribute this dis- 
crepency to a previously unobserved population of z =1-3 Luminous Infra- Red 
Galaxies (LIRGs) with IR luminosities in the range 10^^ Lq < Lir < 10^'^ Lq 
(Shupe et al. (2008); Pearson et al. in prep). 

At longer wavelengths the sensitivity of Spitzer is not great enough to probe 
the high redshift 70^m &: 160/im populations and hence the number counts at 
these wavelengths are generally much closer to those predicted by pie- Spitzer 
models (see Lagache et al. 2004; Frayer et al. 2006; Frayer et al. 2009). 

Finally Le Borgne et al (2009) have used a combination of the published 
number counts at all Spitzer MIPS wavelengths, in combination with the pub- 
lished 15/xm counts from ISO and the 850/xm counts from SCUBA to find the 
best-fit star formation history of the Universe. Their approach involves a non- 
parametric inversion of the counts at the 5 wavelengths. Specifically they blindly 
construct bolometric IR luminosity functions which can match the observed 
number counts, assuming the fliix at each wavelengths is given by an phenomeno- 
logical set of far-IR SEDs. The range of IR luminosity functions and redshifts 
which are allowed by the data are then collapsed to give the range of plausible 
star formation rate densities (SFRD) of the Universe. They find good agreement 
to direct measures of the SFRD as compiled by Hopkins & Beacom (2006) , with 
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the conclusion that LIRGS at z =1-2 and Ultra Luminous Infra-Red Galaxies 
(ULIRGS: < Ljr < W^^Lq) at z ~ 2 are the dominant contributors to 

the SFRD, and hence the GIB. 

While number counts are a useful tool for constraining population models, 
the degeneracy between luminosity and redshift can limit their interpretability. 
Fortunately great effort has been put into large spectroscopic campaigns, as 
well as the production of accurate photometric redshifts, in the deepest Spitzer 
survey fields. 

Taking advantage of data from the VVDS,G00DS,C0MB0-17 and deep 
24/Ltm surveys in the GDFS field, Le Floch et al. (2005) were amongst the first to 
show that the space density of LIRGs increases by a factor of 100 from z =0- 
1. This result has been subsequently confirmed by Magnelli et al. (2009) using 
both MIPS 24/xm and 70/xm data taken as part of the FIDEL survey. At higher 
redshift (z > 2) Caputi ct al. (2007) find that an even greater increase (> 100 
times) in the space density of ULIRGs from z =0-2. These results match up 
well with the picture painted by the analysis of the number counts, with LIRGs 
dominating the far-IR luminosity density, and hence star formation density, at 
z ^ 1 and ULIRGs at z > 2. However it is worth noting that most of this work 
relies strongly on the 24/Ltm data from Spitzer which is probing different parts of 
the mid-IR SED at different redshifts. This is particularly a problem at z > 2 
where the 24/Ltm band is probing ~ 8//m rest-frame, a region which can often be 
contaminated by PAH and molecular line emission. 

3. Properties of star forming galaxies 

While Spitzer MIPS data has been key to uncovering the evolution in star for- 
mation since z ~ 2, by combination with overlapping Spitzer IRAC and ground- 
based optical/near- IR data the nature of galaxies hosting this star formation 
can be determined. 

Using the overlap between the SWIRE survey and the Sloan Digital Sky 
Survey (SDSS), Davoodi et al. (2006) showed the disconnect between the optical 
and mid-far infrared properties of local galaxies (z < 0.2). Specifically they 
found a significant fraction (~ 30%) of red sequence galaxies show an excess in 
their mid infrared emission, defined as log(L24/L3.6) > —0.7). Roughly half of 
this red, mid-IR strong, population are determined to be star-forming systems 
from emission line diagnositics (Baldwin et al. 1981), with the other half coming 
from AGN processes. 

In addition to optical colour (essentially a proxy for stellar age), determina- 
tion of the stellar mass of star forming systems, and its evolution with redshift, 
is of particular interest. 

Zheng et al. (2007) used a combination of COMBO-17 data with Spitzer 
IRAG &: MIPS observations to determine the evolution of the specific star forma- 
tion rate (SSFR: SFR/M©) from z =0-1. They find clear evidence of a significant 
decrease in the SSFR across a wide range of stellar mass (IO^Mq-IO^^Mq) from 
z = 1 to present. While the SSFR is found to decrease with increasing mass, 
as postulated by "downsizing" the evolution of the SSFR is similar across the 
range of stellar mass. Thus it is clear that the mechanism driving this evolution 
is a) different from that at high-z, otherwise massive galaxies would never form. 
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and b) the same for all galaxies at z < 1. Zheng et al. suggest that the decline 
in available infalling cool gas, as opposed to environmental or feedback effects, 
is the underlying cause of this behaviour. 

Damen et al. (2009) and Perez-Gonzalez et al. (2008) have utilised the 
deeper FIDEL MIPS data in CDFS & the extended Groth strip respectively to 
investigate the evolution of the SSFR out to z ~ 2. Damen et al. find good 
agreement with the Zheng et al. results at low-z, again seeing a strong increase 
in the SSFR with redshift across all stelar masses. With the greater depth of 
the FIDEL data they are able to track the SSFR out to z = 2 for the most 
massive galaxies (> lO^^M©). The SSFR in massive galaxies is observed to 
increase by a factor of nearly 10 from z =1-2, qualitatively in line with the 
downsizing paradigm, although again the slope of the evolution in the SSFR is 
found to be independent of stellar mass. Perez-Gonzalez et al. find a similar 
increase in the SSFR with redshift, with an observed increase of nearly 100 x 
from z =0-2. Interestingly they also investigate the morphological dependence 
of this evolution. Seperating spheriods from disk dominated systems via use of 
the Sersic index differential evolution in the SSFR is observed, with the spheriod 
dominated systems experiencing a much stronger increase in SSFR with redshift. 

Using the larger SWIRE dataset, in combination with photometric redshift 
and stellar mass estimates from Rowan- Robinson et al. (2008), Oliver et al. in 
prep investigate the SSFR as not only a function of mass and redshift, but also 
optical spectral type. Figured] shows the results of this analysis. While a similar 
increase of the SSFR with redshift is observed, in contrast with the other results 
some dependence of the evolution with stellar mass is observed, once the sample 
is broken down by spectral type. In particular the elliptical samples appear to 
show a strong decrease in SSFR with increasing mass, while late-type galaxies 
show either no trend or a slight increase. 

4. Environments of star forming galaxies 

It is well accepted that environment plays a key role in the evolution of galaxies 
(e.g. Moore et al. 1996; Dressier et al. 1997; Farrah et al. 2004; Hogg et 
al. 2004; Kauffmann et al. 2004; Bamford et al. 2009.) However actually 
quantifying this role via observations has proven difficult. Spitzer observations 
over wide areas has done much to alleviate this by allowing consistent samples 
of galaxies across a range of environments, redshifts, and spectral types, to be 
assembled. 

Elbaz et al. (2007) investigated the environment of star forming galaxies 
at z ~ 1 in the GOODS field. Using a "counts-in-cells" approach the local 
galaxy number density was estimated by the number of galaxies within boxes 
of comoving size 1.5 Mpc. Combining this information with star formation 
rates estimated from deep MIPS 24/im observations the relationship between 
star formation rate and environment is presented. Curiously a strong increase 
in star formation rate with increasing galaxy number density is observed, in 
stark contrast to what is seen in the local Universe (z ~ 0.1) by surveys such as 
2dFGRS or SDSS. 

Taking a similar, but philosphically different, approach Caputi et al. 2009 
investigated the properties of close neighbours (< 1 Mpc or within 500 kms~^) 
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Figure 1. SSFR as a function of stellar mass and spectral type. Spectral 
types and stellar mass estimates come from the catalogue of Rowan- Robinson 
et al. (2008). 



to 24^m selected LIRGs and ULIRGs at z ~ 0.8 in the COSMOS field. Via 
measurements of the D„(4000) indices of the neighbouring galaxies they find 
that LIRGs lie in regions with a significantly greater fraction of passive galaxies 
than ULIRGs; 42% of LIRG neighbours have D„(4000)> 1.4, compared to 24.5% 
for ULIRGs. Interestingly both of these numbers are significantly different from 
the fraction measured for a control sample of randomly selected galaxies; ~ 30%. 
By looking at the mass-weighted number density around their LIRG and ULIRG 
samples Caputi et al. find that LIRGs at z ~ 0.8 lie in overdense regions, while 
ULIRGs lie in underdense regions. They speculate that this is the underlying 
cause of the difference in the stellar populations of LIRG and ULIRG neighbours. 

To investigate this more quantitatively Frost et al. in prep use the large 
area of the SWIRE survey, in combination with photo-z's and SED fitting from 
Rowan-Robinson et al. (2008), to measure the spatial correlation function for 
SWIRE galaxies as a function of spectral type, redshift and stellar mass. 

Figure [2] shows the clustering strength for early and late-type galaxies as a 
function of redshift and stellar mass. It is clear th a t the colo ur-density relation 
was in place at z = 1 in agreement with Coil et aTl ( 20081 ) and McCracken et al 
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Figure 2. The comoving clustering strength tq as a function of redshift for 
EUiptical (triangles) and Sab-Scd galaxies (squares). Multiple stellar mass 
ranges are shown for each type; IO^^Mq < M, < IO^O-^Mq (light blue), 
IQI'^-^Mq < M, < 10"Mo (green), lO^M© < < 10"-^ Mq (yellow), 
10"-^Mq < < 10"-8Mo (red). 



( 20081 ). Interestingly the clustering of early and la t e-typ e galaxies does not scale 
with stellar mass, in agreement with ICoil et all (|2008|) but at odds w ith low 
redshift measurements iNorberg et al. 1 ([200^) and iBudavari et"all (|2003l 1. This 
implies that at higher redshifts elliptical galaxies are found in the same environ- 
ments regardless of their stellar mass content. Frost et al also find the clustering 
of elliptical galaxies to increase steeply with redshift over 0.1 < z < 1.0 whereas 
it decreases slightly for spirals over the same redshift range. The clustering 
of massive star forming spirals is found to increase strongly with redshift for 
z > 1.0, in qualitative agreement with the work discussed above. 

At high-2 Farrah et al. (2006) use the effect of redshifting the 1.6^m stellar 
"bump" feature through the IRAC bands to isolate samples of galaxies in the 
SWIRE survey at z ~ 1.7 and z ~ 2.5. Farrah et al. select Bump-2 & Bump-3 
sources (sources which "peak" in the IRAC 4.5;um and 5.8/im band, respectively) 
which are also found to be brighter than 400/iJy at 24;um. These selections 
roughly translate to a stellar mass limit of lO^^M© and IR luminosity limit of 
10^^ Mq i.e. the most massive ULIRGs at high-z. Via measurement of their 
spatial correlation, both samples are found to have moderately strong clustering 
strengths, with a clustering length of tq = 6.14 it 0.84 Mpc for Bump-2's 
(1.5 < z < 2.0) and ro = 5.36 ± 1.28 /i^^ Mpc for Bump-3's {2 < z < 3). 
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Interesting this compares well with measurements of clustering of high-z QSO's, 
as shown in Figure [3j Assuming a realistic model for dark matter halo evolution 
to z = shows that these galaxies will occupy high density environments (> 
IO^^Mq dm halo mass). Thus these galaxies are the likely progenitors of > 
ellipticals at z = 0. 
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Figure 3. Comoving correlation length (ro) vs. redshift for the Spitzer 
selected high-z sources as well as a compliation of other measures (Farrah 
et al. 2006 and references therein). Assuming a range of dark matter haloe 
evolution it is clear that high-z ULIRGs evolve to occupy massive (> lO^^Mo) 
DM haloes. 

Similar levels of strong clustering are found for analogous samples of high-z 
ULIRGS. Brodwin et al. (2008) measure the correlation length of Dust Obscured 
Galaxies (DOGs) at 1.5 < z < 2.5 to be tq = TA'^Q^'i. They find the clustering 
strength of DOGs is strongly dependent, with the most luminous sources (F24 > 
O.GmJy) having a correlation length of tq = 12.97I']2'q4. Similarly Magliocchetti 
et al. (2008) find a correlation length of ro = 8.5]'lj^'| for a photo-z selected 
sample of z > 1.6 F24 > 0.4 mJy sources. Thus it is clear that ULIRGs at z ~ 2 
inhabit regions which are the progenitors of z = clusters and groups. 

5. Conclusion 

While the results outlined here form an impressive resume for the achievements 
of Spitzer they represent only a small fraction of the vast accomplishments of 
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the observatory in this field. In particular highlights excluded from this review 
include important work such as that on local galaxies performed by the SINGS 
survey (Kennicutt et al. 2003), and the vast array of discoveries made via mid 
IR spectroscopic observations with the IRS instrument (e.g. Brandl et al. 2006; 
Smith et al. 2006; Yan et al. 2007; Spoon et al. 2008; Pope et al. 2008; Farrah 
et al, 2008; amongt many others). 

While much has already been learnt from the large survey programs con- 
ducted with Spitzer the true legacy nature of these datasets will soon become 
apparent with projects such as Herschel and SCUBA-2 set to re-observe Spitzer 
legacy fields. While these new datasets represent the next frontier in far-IR/sub- 
mm astronomy the existing Spitzer data will be crucial in interpreting these data. 
In this way, combined with a new generation of large Spitzer surveys undertaken 
as part of the "warm" mission, Spitzer will continue to play a key role in the 
study of star forming galaxies and their evolution for years to come. 

Acknowledgments. IR thanks the conference organisors for their kind in- 
vitation and funding support. 
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